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We investigate the geometry of a critical system undergoing a second order thermal phase tran- 
sition. Using a local description for the dynamics characterizing the system at the critical point 
T — Tc, we reveal the formation of clusters with fractal geometry, where the term cluster is used 
to describe regions with a nonvanishing value of the order parameter. We show that, treating the 
cluster as an open subsystem of the entire system, new instanton-like configurations dominate the 
statistical mechanics of the cluster. We study the dependence of the resulting fractal dimension on 
the embedding dimension and the scaling properties (isothermal critical exponent) of the system. 
Taking into account the finite size effects we are able to calculate the size of the critical cluster in 
terms of the total size of the system, the critical temperature and the effective coupling of the long 
wavelength interaction at the critical point. We also show that the size of the cluster has to be 
identified with the correlation length at criticality. Finally, within the framework of the mean field 
approximation, we extend our local considerations to obtain a global description of the system. 



I. INTRODUCTION 



The understanding of the geometry of systems near a second order critical point is the subject of numerous recent 
works 1^. Most of these works consider dynamics in discrete space (lattice) and try to explain the formation of 
clusters with fractal geometry on the embedding lattice in terms of the scaling properties (critical exponents) of the 
system . In a recent work Q we have proposed a mechanism in order to understand the formation of the fractal 
clusters at T = Tc for systems defined in a continuous space. Based on a scale invariant effective action describing 
the dynamics at the critical point we were mainly interested in revealing how this dynamics leads to the formation 
of critical clusters. A general class of saddle-points of the effective action at T = Tc turns out to dominate the 
configurations contributing to the partition function if we consider the statistical mechanics of an open subsystem 
(cluster) of the global critical system. In the present work we present in more detail and completeness the method 
used in Q| to obtain a consistent picture of the local geometry at the transition point for one-dimensional systems 
and then we apply our approach in order to describe critical systems in higher dimensions. We take into account 
also finite-size effects and we discuss the possibility to use different functional realizations for the order parameter 
characterizing the system at the critical point. Based on the local description of the critical system we propose an 
algorithm, using arguments within the framework of the mean field approximation, to construct the global system 
and to determine its scaling properties. 

The starting point in our investigation is the effective action of a thermal system at the critical point T — Tc, 
specified in d dimensions in terms of a macroscopic field 4> (order parameter) as follows: 

TM = giA^ J d'^a:[i(Vrf0)2 + g,A^\A-U\'+'] (1) 

In eq. (^ 51,52 are dimensionless couplings, (p ^ (volume) and the ultraviolet cut-off A of the underlying microscopic 
theory fixes the coarse graining scale Rc ~ A^^ of the effective system. Eq. (|l|) leads to the standard equation of state 
at r = Tc : 

and the index S is identified with the isothermal critical exponent of the system. The action (^), beeing dimensionless, 
implies: /3 = —{d + 2) and 7 = 2d + 2. Introducing also the dimensionless quantities: (p = ! = ^^i, 

rewrite the effective action (|l|) as follows: 

r,[0]=5i / d%^{\/$ f+gS\'^^] (2) 

The scalar quantity describes in general the density of an extensive physical quantity characterizing the phase 
transition (like, for example, magnetization density or particle density). 
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Let us now mention some examples of theories which belong to the class of physical systems described through the 
effective action (||). 

• 0(N) 3-d effective theory: The action, in the large N limit, for spherically symmetric order parameter in 
the internal 0{N) space, is written as |^: 

where X — This action, for = A^'^^ and x = Ax, has the form: 

It belongs to the general class (||) with : gi = , g2 = §(^77^)^, d ^ 3, S = 5. 

• 3-d Ising model: The effective action rc[cr] - which describes effectively the QCD at the critical point (T = Tc) 

- is written as: 

where the macroscopic field a ^ (length)^^. This action, for a = A^^cr and x — xA, has the form: 

Tc[a] ^xj d^f [i(Va)2 + GA^'-^^^+i] 

It belongs to the general class (j^) with gi — X, g2 — GX^^^ , d — 3. We recall that A = 

Throughout this work we use the convention kb = 1 (Boltzmann constant) and the energy is given in inverse length 
units. We will also use, for simplicity, the notation {(j),Xi) instead of {4>,£i) meaning always, unless otherwise stated, 
dimensionless quantities. 

The paper is organized as follows: In Section II we investigate the Statistical Mechanics of the critical system for 
d — I {d is the embedding space dimension) described by an effective action of the form (||). The formation of fractal 
clusters is shown and the corresponding geometrical characteristics (size, dimension) are determined. In Section III 
we extend the analysis to higher dimensions. In Section IV we study the dependence of the geometrical properties of 
the critical clusters on the coarse graining scale of the effective theory. We also apply our approach to critical systems 
with a more general functional form of the order parameter. Taking care of the finite size effects we determine the 
correlation length in terms of the size of the formed clusters. Then using a mean field approach we construct the 
global system as a superposition of individual clusters and we explore its scaling (geometric) properties. Finally in 
Section V we summarize our main results and we give a brief outlook. Some lengthy formulas referred to in the main 
text are given in the Appendix. 



II. THE 1-D SYSTEM 



The statistical mechanics of the critical system is determined through the partition function: 

Zi= J I?[0]e-^^[*1 (3) 

The local description implies that the integration measure in eq. (^ is over field configurations defined in an open ball 
with radius R and center x, subset of the space V (which in fact can be infinite) occupied by the entire system. 
We call a cluster C the set consisting of points belonging to ft {C C fi) for which the order parameter cj) is greater 
than or equal to a minimum value (cut-off) (j>min- We identify then x with the center of the cluster C. Without loss 
of generality we can set a; = 0. The local geometrical properties of the system are determined through the scaling 
properties of the extensive quantities characterizing the cluster C as we vary the radius R. In order to illustrate our 
method we will first consider, for simplicity, the one-dimensional case. The extension however to higher dimensions, 
as we will see in the next section, is straightforward. 
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In the one-dimensional case the effective action of the critical system (T = Tc), resulting from eq.(H), is given by: 

rM^9i£dx[^i^f+g,\cb\'+'] (4) 

Here we will consider models for which the condition gi ^ 1 is valid. This requirement allows us to use the saddle-point 
approximation to evaluate the path integration in eq. by replacing it through a summation over the saddle-points 
of the action (^). As the subsystem f2 is open no boundary conditions restrict the configurations which contribute 
to the path integral (^. This point of view is essential in our approach. 

The saddle-point configurations (f>{x) fulfil the Euler-Lagrange equation — -^^^ where U (0) is the concave 

pontential U{4') — —g2\(f>\^^^- Considering this equation as the motion of a classical particle we get the first order 
equation: 

E=li^)"-9M'^' (5) 

where E' is a conserved (during the classical motion) quantity identified with the total energy of the moving particle. 
Eq.(||) can be integrated (for details see the Appendix) to give, for E — 0, instanton-like solutions of the form: 

0(x) =Ai|x-a;„r^ ; A, ^ {S - ij^ ^ (6) 

with Xo = (</^(0)) • Thus, for E ~ 0, the position of the singularity Xo depends only on the initial 

condition 0(0). For E ^ the solution has the form (j){x) = Ai\x — x'^]^"^ where now x'^ = x'^{(t>{0) , E) (see the 
Appendix). However, configurations with E ^ 0, contribute to the partition function Z with a suppression factor 
g-giR\E\^ suggesting that the dominant saddle points in the path summation (|^) are those solutions of the equation 
of motion for which E Ri 0. In this case eq.(||) simplifies to: 



rc[0] = 25132 / dx i(t){x)) 
Jo 
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Only configurations with Xo > R give a nonvanishing contribution to the path integral (^. In fact, the partition 
function is dominated by those saddle-points for which Xo ^ R- Since x £ (— -R, R) we can easily take: 4>{x) = 

const =^ AiXo It is straightforward to show that these solutions correspond to the long wavelength modes of the 
field (/)(x) by taking the Fourier transform of eq.(^). We get f{k) ~ (^g^^-^j^s-s and the envelope of /(fc) is given by 

Xo 

Using the above approximation we have: 

r, = GiRxo^^ (7) 
with Gi = 2gi52^i^^- The summation over the saddle points of the action (El) becomes, within this approximation, 

_5 + l " 

an ordinary integration over Xo with measure: Dtj) — dfi{xo) ~ Xo dxo- As stated above the singularity Xq must 
be located outside the cluster C to give a nonvanishing contribution to the partition function of the one-dimensional 
system. This condition fixes the lower limit in the integration over Xq to be Xo.miii = R- To determine the upper limit 
of Xo one has to go back to the definition of the cluster C given previously. Without loss of generality we can take the 
extensive quantity characterizing the geometry in C to be the magnetization M = (j>{x)dx fulfilling the condition 
M > fi with /i = (jjjnindx — R(j)min- In the approximation of constant configurations for the order parameter (j) we 
obtain the upper limit for Xo as: Xo < ■ 

The one-dimensional partition function in now written as: 

/•(^)^ _i±i 
Jr 

Using this expression it is straightforward to calculate the mean value of the magnetization in the cluster C: 
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Using (g) we can show analytically (see Section IV) that in the large Gi limit (Gi ^ 1) there are three characteristic 

s+1 1 ^-1 

regions determining the behaviour of the integral in eq.(@). Setting Rd = Gj and i?„ = Gi^^ we find that: 

• For the region i?^ ^ i? ^ i?„ we have: < L (t){x)dx >~ i?^, with coefficient ai w AiGj^ ^rTzilT' leading 
to a fractal structure of the cluster around the point x = with fractal mass dimension [^[jsj: dp = -g^- 

• This behaviour crosses over for R ^ Ru to a different power law: 

< /p (j){x)dx >^ R"^ suggesting the presence of a fractal with mass dimension dp = |5f at large scales. 

• The lower limit Rd defines a minimal scale of the critical system, below which the fractality is broken. 

The fractality in the central region characterizes the critical system in the sense that it corresponds to the scaling 
behaviour in the vicinity of the local observer when fi —^ 0. The crossover scale i?„ gives a measure of the correlation 
length of the finite system aX T = T^. In Fig. la we show the numerical results for the calculation of < M > using 
the values Gi = 5-10^ and (5 = 5. We recognize the central region and the two scales Rd, Ru- In Fig- lb we plot, for 
the same values of Gi and (5, the results for the mean magnetization if we ignore the breaking at Rd (/^ — > 0). 

III. FRACTAL CLUSTERS FOR D > 1 

Let us now extend our approach to higher dimensions starting from the two-dimensional case. We write the efi^ective 
action (||) for d = 2: 

and look for classical saddle-points in an open subset of TZ^. The Euler-Lagrange equation, in this case, has the form: 

V^^^ ((5 + 1)520'' 
and the corresponding instanton-like saddle-points are: 

d=2 , (^2(r) =A2|r-r;r^ ; = {S - 1)^ (S + 1))- ^ (9) 

We proceed in a similar way as for the one-dimensional case considering the partition function Z2 = J I?[(^]e~'"=['^l . In 
the path summation contribute, similarly to the 1-d case, saddle-points for which fo lies outside the cluster C. The 
main effect in the statistical mechanics of the 2-d system is obtained through the summation over paths with \ro\ ^ R 
(i? is, once again, the radius of C) which are in fact constant configurations determined by the 2-d parameter fo- In 
close analogy with the one-dimensional treatment we write the path integral in Z2 as an ordinary integral over fo- In 
this regime (constant configurations) the two-dimensional effective action is : 

F, = G2i?2ro ^ 

with G2 — T^giijj^fyi + .92^2'''^] ~ 27r(/i(72^2~''^ (^4^)- Performing the calculation of the mean value of the mag- 
netization < M{R) >—< J d'^f(j)(f) >, characterizing the two-dimensional critical cluster and using the notation 
i?^ = yl2"^/i^G|' and i?„ = G^ we find: 

• For i?d < i? < Ru- 

< M{R) >- i?^ (10) 
with coefficient a2 ~ 7rA2G, /t^ . This suggests the formation of a geometrical stucture in C with fractal 

^ ( S+l ) 

mass dimension dp — j^- 

• This behaviour crosses over for R^ Ru to a, power law < M[R) >~ R describing a local fractal with mass 
dimension dp — 2|5y at large scales. 
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• Finally, as in the 1-d case, for R <^ Rd the fractality is broken and the mass dimension coincides with the 
embedding dimension. 

The extension to dimensions c? > 3 needs more care. In this case we must take into account the relation between 
the isothermal exponent S and the anomalous dimension r/: S = ^-i+l iMl- ^^^^ examine the case f] = 0. 

Repeating the procedure followed in the 1-d and 2-d case we obtain analytically the saddle-points of the d— dimensional 
critical effective action: 

^ = Adirl - r')"^ ; = (11) 

and transforming the path summation in the partition function Zd into an ordinary integration over ro we find: 

<M{R)>^< jd'^r(t)>^ R^+'^Z'^ (12) 

for Rd R < oo. This means that for d = 3 i?^ ^ oo. Namely, the crossover to the second fractal has disappeared. 
What happens now if we take f] into account? Consider the case < 77 < 1. For a wide range of universality 
classes, including the 0(4) theory where r/ w 0.034 the anomalous dimension i] obeys this condition. Actually 
for 3-d systems rj is very close to zero [|ll| . The corresponding Euler-Lagrange equation: 

V^0= ((5 + 1)^20^^^ (13) 
cannot be solved exactly. Only an approximate instanton-like solution can be obtained analytically: 

Mr) = AdiK ~ r^) ^ ; Ad ^ C ) ' ( ) ' (14) 

Details concerning this calculation are given in the Appendix. Based on the solution (^4|), and following the process 
applied for one and two dimensions, we determine < M{R) > for spherical symmetric clusters in c? > 3 dimensions. 

Usmg Rd ^ a -dr^iirG^' and Ru = + with a = {-^) = , Gd = dT(d/2) 9i92 and q = 2"^, we 

obtain: 

• For i?d < i? < Ru. 

< M{R) >- (15) 

with coefficient ad ~ dV(d/2) ^"f(zil) ■ 

• For R ^ Ru the power-law: 



< M{R) R^ 



d{2-r,) 



Breaking of the fractality for R Rd- 



Comparing eqs.(12) and ([l5| ) we find, for 77 = 0, the same power law. This serves as a consistency check of the 
aproximation we used. We calculated the saddle-points of eq. ( p^ ) numerically. We also have calculated, in the constant 
configuration regime, the mean magnetization < M{R) >. The results are presented in Fig. 2. In Fig. 2a we plot 
together the numerical and the approximate solution to the Euler-Lagrange equation for d = 3 and rj = 0.34 The 
characteristic behaviour of < M{R) > for c? = 3 is presented in Fig. 2b. Here we have used G3 = lO'^ and rj — 0.34. 
The breaking of the fractality (for R <^ Rd) is clearly reproduced while the crossover is suppressed due to the small 
value of 77. 



^We used this value istead of = 0.034 valid for 0(4) in order to magnify the difference between the approximate and the 
numerical solution. 
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Putting together our results for 1, 2 and 3 dimensions and denoting by dp the fractal dimension in the central 
region of the cluster C while by dp the fractal dimension beyond the corresponding upper limit we have found: 

= — — ; d= 1,2,3,.. (16) 



dp^d^- — - ; d=l,2 (17) 

— 1 

While for d > 3, we have: 

dp - dp = ^^^^^ + 0{v') (18) 

A remarkable property of the geometry of the cluster is that dp > dp for all dimensions, indicating a dilution of the 
cluster C for distances greater than Ru from the center of the cluster. In other words, the size Ru of the cluster gives 
a measure of the correlation length in the finite system. For 3-d systems, however, the maximal size of a single cluster 
(Ru) coincides practically with the size of the whole system (77 « 0, dp « dp) and one recovers the conventional 
behaviour of the correlation length ^ in a second-order phase transition (^ ~ size of the system). For critical systems 
of low dimensionality {d < 3) the association of the correlation length with the size of the system needs particular 
care and this issue will be discussed in detail in Section IV. 



IV. EXTENSIONS AND FINITE-SIZE EFFECTS 



In Section III we have shown the appearance of a fractal geometry for the cluster C in the central region of scales 
Rd <^ R Ru- For R < Rd we obtain the breaking of fractality and beyond Ru a more dilute fractal emerges. 
Therefore the limits , Rd determine the part of the cluster with fractal dimension dp . In the following we will 
investigate how a change of the coarse graining scale A affects the fractality region. 

We consider the transformation: 

A-^ = AA'"^ (19) 
where the ultraviolet cut-off A fixes the coarse-graining scale (A^^). Then eq.(||) becomes: 

TM = g,X'+'A'-^^+'^ J d'xi^iVd^f + ff2A-(2rf+2)A''''+'|A'^A'^''0|^+i] (20) 
Setting (j)' = A'~'^(j), x' — A'x eq.(pO|) simplifies to: 

r,[0] ^giA'^+^l d'x'[^{^d4>r + 92X'^'-'^+^'-'^\^'\'+'] 

where the new constants g'l, g'2 have the values: 

9i = ^ 



92 = .92 



^d{S-2) + {d-2) ^21) 



5-1 5- 

We have seen that the dimensionless values of Rd for the 1-d case are: Ru = Gj^ ^ 9T^ 92^ ^"^^ = 

S+l 1 1 1 

~^ Gj ~ 9i 92 ■ According to eq.(21) the new values for the dimensionless limits are: R'^ = XRu and i?^ = XRd- 



Using eq.([l9|) we find that the quantities A ^Ru, A ^Rd do not depend on the choice of the cutoff A: 

A' i?^ =^ A ^Ru 

A'-'r'^ - A-^Rd (22) 

■S-l .S-l _1 s + l 1 1 1 

We may now extend these calculations for the case d = 2 where Ru = G'2 ~ 5i ff2 ^ ^'^'^ = ^2 ~ ffP 9? 



and we get again eqs.(|22|). More complicated is the case d > 3. If we neglect the anomalous dimension 77 we recover 
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eqs.d23) but taking 77 into account we cannot find an analytic expression for the limits Rd, Ru- Using, however, the 



approximate solutions (|14| ) we can prove the validity of eqs.(22) to a leading order in 77. 

Let us now consider the case when the order parameter is not directly the scalar field (j){x) but a power: (/>"(a;), 
n > 0. The extensive variable characterizing the critical geometry is now taken to be: 

M{R) = [ d'^x 
Jc 

Performing the calculation of < M{R) > at the level of the saddle-point approximation we obtain: 

• In the central fractality region R'-f < i? < rI^'^ the dimension is: 

where the embedding dimension d takes the values: d — 1,2, 3, ... 

• The geometry in the external region R ^ R^f' is described through the dimension: 



df^d-^^ ■ d^l,2 (24) 



valid also for d > 3 if we neglect the anomalous dimension rj. 

upper lii 

d{S-l) 



Using the obvious condition dp^ > eq.(^J) leads to an upper limit for the value of n: 



n < 

~ 2 

Using as an example the parameter values d — 1 and 5 = 5 we obtain n = 1, 2 as the possible values of the power 
n. For this special choice of d and 5 we obtain a remarkable property when the order parameter is (j)^: in this case 
dS^p) — and therefore r'^u^ coincides with the correlation length. Concerning the limits r'"^'^ and R^f"^ and their 
dependence on the power n we find: 

• The upper limit R^'^ does not depend on n. 



The lower limit Ru has the following form for a general n: 



We turn now to the question of critical cluster formation in a finite system. The effective action (H) is in fact 
valid for the ideal case of an infinite system. In order to take finite size effects, in a consistent way, into account, we 
have to add the term ^m^cfP in @. In the following we will consider the statistical mechanics, within the theoretical 
framework developed so far, of the modifed effective action, which includes the above mentioned quadratic (mass) 
term in 0. Thereby we restrict ourselves in the simplified 1-d case although our results can be extended to higher 
dimensions without difficulty. We also use for the isothermal critical exponent the value i5 = 5. The central interest 
in our investigations is to determine the changes induced to the upper limit Ru of the central fractality region, due 
to the presence of the mass term. This may lead us to a better understanding of the relation between i?„ and the 
correlation length of the critical system. 

The saddle-points of the modified action are obtained (in the 1-d case) using the energy integral (||) . The dominant 
configurations are those for which E = 0. We have then: 

/■'^'^^^ 







Solving eq.(25) for (j) we get 



Mme^"^- _ 0^(0) 

[m zg2C e ) ' 1 + , /1 + ^04(o) 
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The solution (pq) for small m simplifies to: 



^ix) ^ J - ..r/^ ; x„ = -I''' (27) 

and taking the limit m ^ we recover eq.(||) for (5 = 5. 
The position of the singularity in eq.(|2^) is: 

1 

io = ±— Zn— ^ (28) 

It is easy to show that Xo Xo (see eq.(^) for m ^ 0. In the following and up to the end of this section we will, for 
simplicity, drop the tilde over Xo- Whenever Xq appears in the following formulas it means the expression (p^) . 
Inserting eq.(|2^) in the solution ( p6| ) we finally obtain: 

(j){x) = (W— m)i/2- ^— (29) 

Repeating the arguments of Section II we per form the path summation in the partition function of the finite system 
using the constant solutions, deduced from (P9|), for Xo ^ x: 

0(x) = ( , l^m) (30) 



Within this approximation the effective action of the finite system is: 



e 



r,(i?,Xo) = Gi?^^-— ; G = 2'^'g,g,'/'m^ (31) 

The path summation in the partition function goes over to an ordinary integration over Xo with measure obtained 
from eq.(^): 



Z - / dxo [TO3/2e-"^°(l-e-*'"^°)-i/2-m3/2e-5"^°(l-e-4"^°)-3/2] 



(32) 



In the limit m ^ eq. (|32|) becomes: 



recovering the expression for the partition function Zi found in Section II. 

GR 
Ssinh-^ (2mz) 



Setting now u;(z) = we get from eq.(||): 



dtf{t)e-^ (33) 

(2:0, max) 



where fit) is given by: 



m = mV2[(^)i/a _ (^1 + 1(^)2/3 „ 1 (^)i/3)(^)i/2] 

(£^r^(G^)i/3 (34) 

The value of i?„ is then determined through the condition uj{Ru) = 1. In the limit m —> we find i?„ = 2^''/^gy^(72 
in accordance with the results obtained in Section II. For m 1 we obtain an analytic expression for Ru : 
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„ Inm 

Ru = ^ 35 

For general m however i?u can be determined only numerically. In Fig. 3 we present the results for the dependence 
of Ru on m solving numerically the equation ijj{Ru) = 1. We have used G = 1. The quantity ^ is actually the linear 
size of the critical system. For a large system (m^^ 3> 1), the cluster size i?„ gets saturated, becoming independent of 
m. The long range correlation, in this case, is generated by succesive convolutions of neighbouring clusters. In other 
words, the picture for the global system, emerging from our results (Fig. 3), is a superposition of fractal clusters with 
finite size, which, by coalescense, may create long-range ordering in the critical system. In order to illustrate this 
effect, we have constructed, by simulation, a global system in 2-d as a set of softly interacting clusters with prescribed 
geometrical properties. In fact, based on a mean field approximation scheme, it is straightforward to determine the 
distribution of N clusters, P{Ri, ...Rn), with centers located at the points {Ri...Rn) in the area Sg of the global 
system. For this purpose we consider the pontential term in the effective action: 

C/(0) = 3152 / '[0(r ; ^^^< f d^f '</>(f ') > (36) 

where Sr is the area occupied by a cluster of radius r and the mean field ^(f) is written, according to our results in 
Section III, as follows: 

r(— ) 

m = ^7rrTi^9i92i6 + 3)/2)-^r-7TT (R^ < r < Ru) (37) 



From eqs.(36) and ( p7| ) we obtain: 



27r 

(38) 



where R{d) specifies the distortion of the area occupied by the cluster in question, owing to a state of coalescense 
with neighbouring clusters. Introducing the mean radius R = R{9) in the integral (p8|) we finally obtain: 



UW = U^^] fT^Vnf^) ; Sd^nRl ; S = nR^ (39) 



where S* is a measure of the area occupied by a distorted (in general) cluster. The distribution of A'^ clusters, in this 
picture, P{Ri...Rn) ^ Y[i=i e~'^*('^\ is given by the following area-law: 



qNas p 

Zn = d^Ri...d^RNiSiS2...SNy 



The smallness of the exponent as in eqs.(pO|) {as ~ 0.006 for 5 — 5) guarantees that the interaction of clusters is 
very weak, leading to a random distribution P{Ri...Rn) in the area Sg and to a Poisson behaviour of the partition 

function, Zjv ^ ^d"'' ih- permits us, therefore, to treat the global system as an almost ideal gas of clusters. The 
parameters of the effective theory as well as the size of the critical system, determine then the number of formed 
clusters, the density within each cluster and the size of each cluster. 

It is now straightforward to construct the global system as a set of softly interacting clusters with prescribed 
geometrical properties. For simplicity let us consider a system for which the order parameter is interpreted as density 
of particles. Given the linear size R, the isothermal critical exponent S and the effective couplings gi, 52 we calculate 

the size i?„ of a single cluster and the corresponding number of clusters Nd — (iT") system. The number of 

particles n within each cluster is then given through eg. (pO|) . The entire net of the clusters is constructed in two steps. 
First we generate the positions of the centers of the clusters treating them as uniformly distributed random variables 
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over a square with site R. Then we generate the points inside each cluster with a distribution function specified by 
eq. (|l^) . If two clusters (say the j-th and the j-th cluster) overlap then a point in the i-th cluster is taken into account 
if the nearest center to it is the center of the i-th cluster otherwise this point is neglected. Such a construction is 
presented in Fig. 4. We have used the parameter values: R = 1, S = 5, gi = 50 and 32 = 1- 

Although the fractal mass dimension of each cluster is the same the resulting global set turns out to have a different 
fractal dimension. In fact calculating the generalized dimensions of the global set we observe that its geometric 
structure does not cor resp ond to a pure monofractal set. Here we have used the method of factorial moments to 
perform this analysis |1J]. We divide the region of the global system into boxes of linear size I {M = y). 
Denoting by rii the number of points within the i-th box we define the p-th order factorial moment of the distribution 
of the points of the global set in space as: 

■j^'^ni{ni - l)..{ni - p + 1) 
(Af ) = (41) 

For a fixed value of p the moment F^p\M) possesses a power-law dependence on M (for M ^ 1) provided that 
the point-set under consideration has a fractal structure: F^p\M) ^ APp . The exponent Sp is related to the 
fractal dimension of the corresponding point-set: Sp = (j) ~ l)((i — dp). The dimensions dp are the generalized 
dimensions characterizing the system d2 beeing the corresponding average fractal dimension. For a monofractal set 
the generalized dimension spectrum is given as: dp — c?2 for p = 3,4, ... We have calculated the first three factorial 
moments {p = 2, 3, 4) as a function of M for the set shown in Fig. 4. The results are shown in a log-log plot in 
Fig. 5. We find the exponents (slopes in the log-log plot): S2 = 0.65, S3 = 1.61 and 54 = 2.6 suggesting that the 
underlying set is a multifractal. A deeper understanding of the dimension spectrum of the global system based on 
the construction described above is a subject for a future investigation. 



V. CONCLUSIONS 



We have studied in detail the formation of critical clusters in a wide class of systems, undergoing a thermal phase 
transition of second order. Wc have used suitable, instanton-like, saddle-point configurations for the local field- 
fluctuations, in order to saturate the path summation of the partition function. In this treatment we were able to 
describe the critical system locally and specify the geometrical properties of a single critical cluster in terms of the 
parameters of the effective action, at the critical point. Our main results are summarized as follows: 

1. In critical systems of low dimensionality {d < 3) there exist two characteristic scales {Rd, Ru) which specify the 
geometry of any single critical cluster, namely its maximal size {Ru) and its fractality region Rd <^ R ^ Ru- 
The fractal dimension is dp = in agreement with other treatments on a lattice, and the minimal scale i?^, 
below which fractality breaks down, is the analogue of the lattice spacing in any treatment in discrete space. 

2. Beyond the scale Ru {R ^ Ru) the fractal dimension crosses over to smaller values, dp < dp, and for a suitable 
choice of the order parameter {(j)^{x) with n = '^^^^^'> ) it may even vanish {dp = 0). Therefore the scale Ru 
can be associated with the direct correlation length, sa RuA'~^ |9|, which coincides with the maximal size of 
a single cluster. 

3. For sufficiently large systems (size » A^^) the direct correlation length {^d ~ RuA~^) is independent of the size 
of the system (Fig. 3) and remains finite even in the thermodynamic limit (infinite system). 

4. The global system is built-up by a random distribution of critical clusters which may overlap giving rise to 
a long-range total correlation and therefore to density fluctuations at all scales. We have shown the validity 
of such a mechanism in 2-d by developing a suitable algorithm in order to construct the global system, based 
on the local description. Our results show that the entire critical system develops strong density fluctuations, 
of multifractal nature, in a wide range of scales, far beyond the size of a single cluster (Figs. 4,5). A deeper 
understanding of this global structure and in particular of the connection between the fractal geometry of a 
single cluster and the multifractal spectrum of the entire system, remains a challenging, open question. 
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5. In 3-d critical systems the anomalous dimension is approximately zero (ry « 0) and the crossover scale Ru 
associated with a single cluster coincides practically with the size of the global system {dp ~ dp)- This 
observation leads to a simple picture according to which the development of fluctuations at all scales, at the 
critical point, is realized through the formation of self-similar clusters at all sizes. The maximum cluster size 
{Ru) coincides with the size of the global system and gives a measure of the correlation length which becomes 
infinitely large, in the thermodynamic limit (infinite system). In this case the geometry of the global system 
coincides with the geometry of a single cluster and therefore it remains monofractal with the same fractal 
dimension dp. 

Concluding, we have investigated the geometrical structure of critical fluctuations, developed locally in a thermal 
system which is described effectively by a scalar field. The formation of fractal clusters with a mass dimension dp, 
fixed by the universality class, has been revealed and the mechanism for generating fluctuations at all scales, in 
the entire system, based on the local properties, has been discussed. It is of interest to note that the fluctuations 
of the global system obey a different geometry in 3-d critical systems (monofractal) and in critical systems of low 
dimensionality, d < 3 ( mult ifr act al ) . The reason for this different behaviour is due essentially to the fact that in 3-d 
systems, described by scalar theories, the anomalous dimension rj turns out to be very small (ry » 0). 

APPENDIX 



Integrating the first order differential equation (||) we obtain: 



d<i) 



^7^(0) [2(£;-f 320^+1)11/2 



It follows that: 



where: 




(E. 1) 



a = 5 + 1 



92 



l/a 

( 

l/a 



^2= (I) %(x) (E.2) 
The integral / can be determined analytically: 



(i + ^)^'^"io (TTvO^'^ 

= a[r2 2Fi(l/2,l/a;l+i,-r^)-ri 2^1(1/2, l/a; 1 + i, -r?)] (E. 3) 



where 2^"i is the hypergeometric function. 
Inserting the formula p3] : 



2F,{a, P; 7; z) = "i (-l)"z-" 2Fi{a,a + l-r,a + l-P; -) 

r(/3)r(7-a) z 



in (E. 3) we obtain: 
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Cl 



(E. 4) 



with: 



and 



Cl 



r(i/a)r( 



2^ a' 



5-1 



f, = ^F,{ll2,\--,\~-,--cj>{xr^) 
z a Z a g2 

/2 = 2^^i(l/2 i--,|-i,--<?i(0)-«) 
2 a 2 a g2 



Then ( [E. 4) ), for 0, leads to the sohition (|6|): 



(-5 - 1)^32 



with 



Xo = Xoi4>iO)) 



-.{my 



From ( [E. 4| ) and for 7^ we obtain that: 

"(^-1)^52 



{x',-x) ^ 



fl 



where now xj, = Xof2{E). Namely, if i? 7^ 0, Xo depends on two parameters: 

x',^x',{m,E) 



(E. 5) 



(E. 



(E. 7) 



(E. 8) 



(E. 9) 
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Figure Captions 



Figure 1: (a) The mean magnetization < M{R) > as a fmiction of B for a 1-d critical system. The parameters are 
chosen so that: Gi = 5 • 10^ and fj, = 1. A hnear fit is also shown in order to indicate the two different fractality 
regions described in Section IL (b) The same plot as in (a) but now with /i ~ 0. The scale Rj, for the breaking of the 
fractality is in this case absent. All the presented quantities are in arbitrary units. 

Figure 2: (a) A typical 3-d instanton-like saddle-point for 77 = 0.34. Both the analytic approximation (dotted line) 
as well as the numerical calculation (solid line) are shown, (b) The mean magnetization < M{R) > for the 3-d case 
using the saddle-points of the form presented in (a) to perform the corresponding statistical averaging. As in Fig. 1 
the displayed quantities are in arbitrary units. 

Figure 3: The upper limit i?„ of the central fractality region for (i = 1 as a function of the size ^ of the critical 
system. We use arbitrary units for R^ and m. 

Figure 4: The global 2-d critical system described through an effective action of the form (j^). We used the parameters: 
(71 50, 52 = 1 and (5 = 5. Each full circle represents a point x in the 2-d space occupied by the critical system with 
<j){x) > (j)min- The coordinates x and y are given in arbitrary units. 

Figure 5: The first three factorial moments for the point-set presented in Fig. 4. The linear fits indicate the slopes 
of the corresponding moments suggesting the multifractality of the underlying set. 
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